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LONG-TERM  GOALS 

The  primary  goal  is  to  demonstrate  the  utility  of  the  ROMS  framework  in  a  real-time,  sea-going 
environment  for  prediction  studies  in  the  Intra- Americas  Sea  (IAS)  with  particular  emphasis  in  the 
Caribbean  Sea.  As  such,  we  will  demonstrate,  as  a  proof  of  concept,  the  utility  of  adjoint  modeling 
and  4-Dimensional  Variational  (4DVar)  data  assimilation  in  a  real-time  operational  setting,  at  sea. 

OBJECTIVES 

The  main  objectives  are: 

1.  To  develop  a  real-time  data  assimilation  and  prediction  system  for  the  IAS  based  on  a 
continuous  upper  ocean  monitoring  system. 

2.  To  demonstrate  the  utility  of  variational  data  assimilation  in  a  real-time,  sea-going 
environment. 

3.  To  to  demonstrate  the  value  of  collecting  routine  ocean  observations  from  specially  equipped 
ocean  vessels. 

4.  To  develop  much  needed  experience  in  both  the  assimilation  of  disparate  ocean  data  and  ocean 
prediction  in  regional  ocean  models. 
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APPROACH 


The  Regional  Ocean  Modeling  System  (ROMS)  will  be  used  for  data  assimilation  and  prediction 
studies  in  the  IAS.  Through  an  innovative  partnership  between  the  University  of  Miami  Rosenstiel 
School  of  Marine  and  Atmospheric  Science  (RSMAS),  the  Royal  Caribbean  Cruise  Line  (RCCL), 
NOAA  and  NSF,  a  RCCL  cruise  ship,  the  Explorer  of  the  Seas,  has  been  equipped  with  a 
comprehensive  suite  of  oceanic  and  atmospheric  sensors  which  provides  continuous  observations 
along  two  cruise  tracks  that  circumnavigate  the  Caribbean  Sea  once  every  two  weeks.  This  is  an 
unprecedented  research  data  set,  and  provides  a  unique  opportunity  to  test  ROMS  and  its  associated 
ocean  data  assimilation  and  prediction  systems  in  real-time  and  at  sea. 

The  IAS  provides  a  non-trivial  test  bed,  and  an  opportunity  to  address  outstanding  physical 
oceanographic  issues  regarding  transport  pathways,  and  mesoscale  variability  in  a  vital  area  of  the 
world  ocean.  The  IAS  is  the  source  of  the  Gulf  Stream  system  and  moisture  from  the  IAS  is  the 
primary  source  of  precipitation  over  much  of  the  eastern  and  central  continental  United  States.  In 
addition,  the  IAS  circulation  exerts  a  strong  influence  on  the  formation  and  propagation  of  tropical 
depressions  and  hurricanes  that  threaten  the  United  States  and  Caribbean  nations. 

The  ocean  state  estimation  and  data  assimilation  methods  used  in  this  project  are  based  on  4DVar.  A 
suite  of  4DVar  data  assimilation  platforms  will  be  used  in  this  project  including  strong  constraint, 
incremental  4DVar  (IS4DVAR),  an  indirect  representer-based  weak  constraint  4DVar  (W4DVAR), 
and  a  weak  constraint  Physical  Space  Analysis  System  (W4DPSAS).  W4DVAR  is  based  on  the 
Oregon  State  University  Inverse  Ocean  Model  (IOM)  of  which  ROMS  is  also  a  component  via  a 
separate  NSF  funded  effort  (Di  Lorenzo  et  al,  2006;  Muccino  et  al,  2006).  The  strong  constraint  4DVar 
assumes  errors  only  in  the  observations  and  initial  conditions  whereas  weak  constraint  4DVar 
considers  errors  in  both  observations  and  model  dynamics. 

WORK  COMPLETED 

During  our  first  IAS  project  PI  meeting  in  Boulder  in  March  2005,  the  following  general  tasks  were 
identified  as  important  steps  toward  achieving  the  goals  of  this  component  of  the  project:  (1)  data 
assembly,  (2)  model  configuration  and  testing,  (3)  experimental  design,  and  (4)  forecast  experiment 
preparation  and  logistics.  These  goals  were  refined  and  reiterated  during  a  second  PI  meeting  in 
Boulder  in  October  2005.  The  following  is  an  up-to-date  status  report  with  respect  to  the 
aforementioned  tasks: 

(1)  Data  Assembly 

(i)  Sub-daily  (swath-by-swath)  satellite  scatterometer  wind  data  has  been  processed, 
assembled  and  used  for  IAS  ROMS  simulations. 

(ii)  NCEP  reanalysis  surface  data  have  been  prepared  and  processed  for  the  period  1990- 
2006  and  used  for  IAS  ROMS  simulations. 

(iii)  Gridded  SSH  and  SST  data  have  been  assembled  and  processed  for  the  IAS  for 
assimilation  into  IAS  ROMS.  The  SSH  data  are  from  AVISO  to  which  a  steric  height 
correction  is  applied,  and  SST  is  a  blended  product  from  several  satellite  platforms 
provided  to  us  by  Dr.  Dave  Foley  at  NOAA  PFEL. 


(2)  Model  Configuration  and  Testing 

(i)  ROMS  has  been  configured  for  the  IAS  region  at  both  20km  and  40km  resolution, 
using  climatological  open  boundary  conditions  derived  from  a  long  run  of  ROMS  for  the 
entire  North  Atlantic.  The  40km  resolution  model  is  designed  as  a  training  tool  for  the 
IS4DVAR  data  assimilation,  and  to  provide  relatively  low  cost  first-guess  initial 
conditions  for  the  20km  data  assimilation  model  runs  and  predictions. 

(ii)  Forced  simulations  spanning  the  period  1990-2006  have  been  run  at  both  20km  and 
40km  and  are  used  as  a  baseline  for  data  assimilation  experiments. 

(iii)  A  series  of  adjoint  sensitivity  calculations  have  been  performed  to  explore  the 
sensitivity  of  the  IAS  circulation  to  variations  in  different  aspects  of  the  model 
configuration.  These  are  described  in  Chhak  et  al.  (2006)  and  Chhak  (2006). 

(iv)  A  series  of  identical  twin  IS4DVAR  experiments  have  been  performed  using  IAS 
ROMS  which  have  greatly  aided  in  the  experimental  design,  detection  of  bugs  in  the 
algorithms,  and  tuning  of  algorithm  performance  and  efficiency. 

(3)  Experimental  Design 

(i)  Semi-variagram  techniques  have  been  employed  to  estimate  the  spatial  correlations 
which  are  critical  paremeters  in  the  background  error  covariance  function  used  in  the 
definition  of  the  IS4DVAR  cost  function. 

(ii)  Based  on  all  of  our  preliminary  experiments,  many  algorithmic  improvements  have 
been  made  to  the  ROMS  4DVar  platforms,  which  include  full  parallelization  for  multi¬ 
processor  computers. 

(4)  Forecast  Preparation  and  Logistics 

(i)  The  proof-of-concept  sea-trials  of  the  IAS  ROMS  data  assimilation  and  ensemble 
prediction  system  will  take  place  aboard  the  Explorer  of  the  Seas,  a  fully  instrumented 
Royal  Caribbean  Cruise  Line  (RCCL)  cruiseship  based  in  Port  Miami.  In  readiness  for  the 
sea-trials,  the  satellite  SST  and  SSH  data  of  l(iii)  for  the  period  March  2005-June  2006 
have  been  assimilated  into  IAS  ROMS  at  14  day  intervals  using  IS4DVAR.  This  will 
continue  through  to  the  end  of  February  2007  when  the  sea-trials  conclude. 

(ii)  A  prototype  real-time  ensemble  prediction  system  is  currently  under  development  and 
will  be  deployed  aboard  Explorer.  Each  ensemble  member  will  have  the  same  initial 
condition  but  will  be  forced  with  different  surface  fluxes  of  heat  and  momentum  generated 
by  a  Bayesian  Hierarchical  Model  (BHM)  that  has  been  developed  by  Milliff  and 
collaborators  for  the  Meditteranean  under  separate  ONR  support.  The  BHM  will  be 
reconfigured  for  the  IAS. 

(iii)  The  real-time  sea-trials  of  the  ensemble  prediction  system  aboard  Explorer  are 
planned  for  the  period  1/7/07  to  2/18/07.  During  this  6  week  period  a  different  project 
team  member  will  travel  aboard  the  cruiseship  running  a  suite  of  ensemble  forecasts  (~100 
per  cruise).  The  initial  conditions  for  each  ensemble  forecast  will  be  generated  using 
satellite  SST,  SSH,  and  IS4DVAR,  and  if  logistics  allow,  using  also  ADCP  and  surface 
SST  and  SSS  data  collected  aboard  the  cruiseship.  As  the  ensemble  unfolds,  each  member 
of  the  ensemble  and  the  ensemble  mean  will  be  compared  against  new  observations  both 
from  satellite  and  the  cruiseship  as  they  become  available. 

(iv)  A  project  website  has  been  constructed  (http://marine.rutgers.edu/po/ias)  which 
describes  the  project  and  acts  as  a  repository  for  the  data  assimilation  products.  During  the 
sea-trials  this  will  also  act  as  a  live  update  for  the  ensemble  prediction  results. 


RESULTS 


Two  different  integrations  of  IAS  ROMS  will  be  compared  here.  In  the  first  case,  referred  to  as  IAS-F, 
the  model  was  run  for  the  period  Jan  1990- July  2006  using  surface  forcing  fields  derived  from  NCEP 
reanalysis  data  and  satellite  scatterometer  wind  data  when  available.  In  the  second  case,  referred  to  as 
IAS-H,  SST  and  SSH  data  were  assimilated  into  the  model  every  2  weeks  using  a  14  day  assimilation 
window  starting  in  March  2005  using  IS4DVAR.  The  first-guess  used  on  1  March  2005  was  the 
model  state  from  IAS-F  on  the  same  date.  At  the  end  of  each  2  week  period,  a  14  day  hindcast  was  run 
using  the  best  estimate  of  the  ocean  state  derived  from  the  IS4DVAR.  All  cases  were  run  using  40km 
horizontal  resolution. 

Figures  1  and  2  show  two  measures  of  the  difference  between  the  model  SST  and  SSH  and 
observations  from  IAS-F  and  IAS-H  for  the  period  March  2005  -  June  2006.  The  measures  used  are 
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Figure  1:  MSE  and  CC  for  SST  and  SSH  computed  within  the 
Caribbean  Sea  for  IAS-F,  IAS-H  and  IAS-A. 

mean  squared  differences  between  the  model  and  observations  (denoted  MSE),  and  the  pattern 
correlations  between  the  model  fields  and  observed  fields  (denoted  CC),  each  computed  separately  for 
the  Caribbean  Sea  (Fig.  1)  and  Gulf  of  Mexico  (Fig.  2).  In  addition  to  IAS-F  and  IAS-H,  the 


corresponding  differences  between  the  model  and  the  observations  during  each  14  day  assimilation 
cycle  are  also  shown  (denoted  IAS-A).  Comparing  IAS-F  and  IAS-H  reveals  that  in  general  data 
assimilation  has  a  positive  impact  on  the  model  hindcasts,  although  this  is  not  always  the  case.  The 
MSE  in  the  Gulf  of  Mexico  (Fig.  2)  for  IAS-H  and  IAS-A  increases  considerably  during  the  period 
April  to  August  2005,  and  again  May  2006  and  exceeds  that  of  IAS-F.  This  appears  to  be  due  to  a 
combination  of  forcing  error  and  model  error.  More  detailed  analyses  reveal  that  the  periods 
characterized  by  strong  hurricanes  (in  this  case  Katrina,  Wilma  and  Rita)  are  also  problematic  in  the 
Gulf  because  the  physical  parameterizations  in  ROMS  are  not  designed  to  handle  such  extreme  surface 
conditions. 

Comparing  IAS-H  and  IAS-A,  Figs.  1  and  2  show  that  the  agreement  between  the  model  and  the 
observations  improves  during  data  assimilation  cycle,  indicating  that  IS4DVAR  is  performing  as  it 
should.  Figures  1  and  2  also  reveal  that  in  general  the  model  and  assimilation  perform  better  in  the 
Caribbean  Sea  than  in  the  Gulf  of  Mexico.  This  is  to  be  expected  since  the  40km  resolution  used  in  this 
case  poorly  resolves  the  Loop  Current  and  its  associated  dynamics. 
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Figure  2:  MSE  and  CC  for  SST  and  SSH  computed  within  the 
Gulf  of  Mexico  for  IAS-F \  IAS-H  and  IAS-A. 


The  pattern  correlations  CC  of  Figs.  1  and  2  reveal  that  despite  the  presence  of  persistent  MSE  during 
IAS-H  and  IAS-A,  the  SST  and  SSH  patterns  agree  very  well  for  much  of  the  period  with  those 
observed.  Exceptions  are  the  Spring- to- Summer  transition  period  in  the  Gulf  of  Mexico  identified 
earlier  (Fig.  2)  when  the  model  has  problems  recovering  the  circulation.  The  very  low  values  of  CC  at 
all  time  for  IAS-F  indicates  that  the  model  forced  by  only  the  observed  forcing  is  unable  to  capture  the 
circulation  patterns. 

An  example  of  differences  in  SST  and  SSH  between  the  model  and  observations  for  one  particular  day, 
29  June  2005  for  IAS-F,  IAS-H  and  IAS-A  is  shown  in  Fig.  3.  Figure  3  reveals  that  considerable 
differences  can  exist  between  the  model  and  observations  in  IAS-F  when  there  is  no  data  assimilation. 
These  differences  are  reduced  dramatically  by  IS4DVAR  (IAS-A),  however  when  data  assimilation 
ceases  and  the  model  is  run  in  hindcast  mode  using  the  observed  surface  forcing  (i.e.  IAS-F),  the 
model  minus  observation  differences  can  increase  considerably  due  to  uncertainties  in  the  surface 
forcing,  model  error,  and  errors  in  the  unobserved  components  of  the  circulation  that  are  not 
completely  constrained  by  the  surface  observations. 
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Figure  3:  SST  differences  (top  row)  and  SSH  differences  (bottom  row)  between  the  model  and 
observations  for  IAS-F  (left  hand  panels),  IAS-H  (center  panels)  and  IAS-A(right  hand  panels). 


IMPACT/APPLICATIONS 


This  project  will  demonstrate  the  utility  of  ROMS  adjoint  modeling  and  variational  data  assimilation 
algorithms  in  a  real-time,  sea-going  environment,  and  will  very  likely  strengthen  the  case  for  installing 
oceanographic  instruments  and  ocean  monitoring  laboratories  onboard  other  commercial  or  Navy 
cruise  ships  plying  different  routes. 

TRANSITIONS 

The  various  ROMS  toolkits  and  4DVar  data  assimilation  algorithms  are  available  from  the  Rutgers 
website  and  are  being  actively  used  and  further  developed  by  other  research  groups  in  the  U.S., 
Australia  and  Chile.  To  facilitate  this  transition  a  mini-workshop  on  the  ROMS  data  assimilation 
platforms  was  held  at  the  University  of  California  Santa  Cruz  in  Spring  2006  and  attended  by  students, 
post-docs  and  other  researchers  from  Rutgers  University,  Scripps,  and  UC  Santa  Cruz 

RELATED  PROJECTS 

The  work  described  here  is  intimately  related  to  the  following  ONR  supported  projects: 

“ROMS/TOMS  Tangent  Linear  and  Adjoint  Models:  Data  Assimilation  Tools  and  Techniques”,  PI 
Heman  Arango,  grant  number  N00014-05-1-0366. 

“Application  of  the  ROMS/TOMS  Tangent  Linear  and  Adjoint  Models  to  the  Littoral  Ocean  and 
Semi-Enclosed  Seas”,  PI  Andrew  Moore,  grant  number  N0001 4-01-1-0209. 

“ROMS/TOMS  Tangent  Linear  and  Adjoint  Models:  Data  Assimilation  Tools  and  Techniques  ”,  Pis 
Arthur  Miller  and  Bruce  Comuelle,  grant  number  N00014-99-1-0045. 

“ROMS  Data  Assimilation  Tools  and  Techniques”,  PI  Emanuele  Di  Lorenzo,  grant  number  N00014- 
05-1-0365. 

“Ocean  State  Estimation  and  Prediction  of  the  Intra-Americas  Sea”,  Pis:  Andrew  Moore  and  Ralph 
Milliff  grant  numbers:  N00014-05-M-0277  and  N00014-05-M-0081 
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